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f**^ ' Abstract. Cold three-flavor quark matter at large (but not asymptotically large) densities may exist as 

^ a crystalline color superconductor. We explore this possibility by calculating the gap parameter A and 

, free energy f2(A) for possible crystal structures within a Ginzburg-Landau approximation, evaluating f2(A) 

K ■ to order A®. We develop a qualitative understanding of what makes a crystal structure stable, and find 

^ ' two structures with particularly large values of A and the condensation energy, within a factor of two of 

' those for the CFL phase known to characterize QCD at asymptotically large densities. The robustness of 

, these phases results in their being favored over wide ranges of density and though it also implies that the 

, Ginzburg-Landau approximation is not quantitatively reliable, previous work suggests that it can be trusted 

I . for qualitative comparisons between crystal structures. We close with a look ahead at the calculations that 

^ ' remain to be done in order to make contact with observed pulsar glitches and neutron star cooling. 

\o : 

, We explore properties of quark matter at zero temperature and at densities that may occur within 

\^ ' neutron star cores. Absent interactions, each flavor of quark will fill momentum eigenstates up to a 

, Fermi surface. However, this configuration is unstable to the formation of Cooper pairs in the presence 

^ ■ of any attractive interaction between quarks. In QCD, the interaction between pairs of quarks that are 

, antisymmetric in color is attractive, and hence the ground state features a diquark condensate that is 

f| ' dominantly antisymmetric in color. This is color superconductivity. We consider condensates antisymmetric 

Oh, in Dirac indices and, consequently, in flavor also, implying that quarks in a pair have different flavors. 

pi • At asymptotic densities {Mg/n 0), where the three light quarks can be treated as massless, quark 

QJ , matter exists in the CFL phase ^ in which quarks of all three colors and all three flavors form Cooper pairs 

' with zero total momentum, yielding (ds), (us) and (ud) condensates, and in which all fermionic excitations 

^ , are gapped, with a gap parameter Aq ~ 10 — 100 MeV. However, at densities relevant for neutron star 

phenomenology, meaning quark chemical potentials at most fi ^ 500 MeV, the strange quark mass Ms 
^ . , cannot be neglected. In neutral unpaired quark matter in weak equilibrium, Ms introduces splitting between 

■ the Fermi surfaces for quarks of different flavor which can be taken into account to lowest order in / fi^ 

by treating the quarks as if they were massless but with chemical potential splittings 6fi2 = {fJ^u — Ms)/2 and 
(5/^3 = (/id — /iu)/2 given by S^2 ~ = = ^'^s/i^^^)■ Note that the splitting between unpaired Fermi 
surfaces increases with decreasing density. In the CFL phase, the Fermi momenta are not given by these 
optimal values for unpaired quark matter; instead, the system pays a free energy price oc Sfji^iJ? to equalize 
all Fermi momenta and gains a pairing energy benefit oc Ag/i^. As a function of decreasing density, there 
comes a point (at which (S/i sa Ao/4 j2j) when the system can lower its energy by breaking pairs. Below this 
density the CFL phase can certainly not be the ground state of matter. Within a spatially homogeneous 
pairing ansatz, the phase that results when CFL Cooper pairs start to break is the gapless CFL (gCFL) 
phase of Ref. [2- However, this phase turns out to be unstable to the formation of counter-propagating 
currents, presumably leading to inhomogeneity, and therefore cannot be the ground state of matter at any 
density j^j. Within the range of densities in which the gCFL phase has lower free energy than the CFL and 
unpaired phases, the true ground state of dense matter must have lower free energy still. We have recently 
proposed a candidate phase (actually, two candidate phases) for the ground state of matter over much of this 
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density regime and at still lower densities 0]. In these proceedings we describe this proposal. For further 
references, in particular to other suggested resolutions of the gCFL instability, see Ref. 0]. 

Crystalline color superconductivity |31 E] naturally permits pairing between quarks living at split Fermi 
surfaces by allowing Cooper pairs with nonzero net momentum. In three-flavor quark matter, this allows 
pairing to occur even with the Fermi surfaces split in the free-energetically optimal way as in the absence 
of pairing [71 151 0| . This is the origin of the advantage that crystalline color superconducting phases have 
over the CFL and gCFL phases at large values of the splitting 6fi. For example, by allowing u quarks with 
momentum p + to pair with d quarks with momentum — p + qg, for any p, we can pair u and d quarks 
along rings on their respective Fermi surfaces In coordinate space, this corresponds to a condensate 

of the form (ud) ^ A3 exp (2iqs ' r) ■ The net free energy gained due to pairing is then a balance between 
increasing jqaj yielding pairing on larger rings while exacting a greater kinetic energy cost. The optimum 
choice turns out to be jqaj = r]Sfi3 with 77 = 1.1997, corresponding to pairing rings on the Fermi surfaces 
with opening angle 67.1° [S]. It is possible to cover larger areas of the Fermi surfaces by allowing Cooper 
pairs with the same jqaj but various qg. yielding (ud) ^ A3 exp(2 i qg • r) with the chosen from 
some specified set {qg}. This is a condensate modulated in position space in some crystalline pattern, with 
the crystal structure defined by {qg}. In this two-flavor context, a Ginzburg-Landau analysis reveals that 
the best {qg} contains eight vectors pointing at the corners of a cube, say in the (±1,±1,±1) directions in 
momentum space, yielding a face-centered cubic structure in position space 

We use the following ansatz for the three- flavor crystalline color superconducting condensate 0: 

A,^exp(2iq?.r) . (1) 

This is antisymmetric in color (a,/?), spin, and flavor indices and is thus a generalization of the CFL 
condensate to crystalline color superconductivity. We set Ai = 0, neglecting {ds) pairing because the d 
and s Fermi surfaces are twice as far apart from each other as each is from the intervening u Fermi surface. 
Hence, I can be taken to run over 2 and 3 only, {qj} and {qg} deflne the crystal structures of the (us) 
and (ud) condensates respectively. We only consider crystal structures in which all the vectors in {qj} are 
equivalent to each other, and same for {qg}, as this justifies our simplifying assumption that the (us) and 
(ud) condensates are each specified by a single gap parameter (A2 and A3 respectively), avoidin having to 
introduce one gap parameter per q. We furthermore only consider crystal structures which are exchange 
symmetric, meaning that {q2} and {qg} can be exchanged by some combination of rigid rotations and 
reflections applied simultaneously to all the vectors in both sets. This simplification, together with 6fi2 = Sfi3 
(an approximation corrected only at order Mf/fi^), guarantees that we find solutions with A2 = A3. 

We analyze and compare candidate crystal structures by evaluating the free energy ri(A2, A3) for each 
crystal structure in a Ginzburg-Landau expansion in powers of the A's. This approximation is controlled if 
A2, A3 ^ Ao, S^, with Ao the gap parameter in the CFL phase at M"^ / fi = 0. The terms in the Ginzburg- 
Landau expansion must respect the global U{\) symmetry for each flavor, meaning that each A/ can only 
appear in the combination |A/p. (The U(\) symmetries are spontaneously broken by the condensate, but 
not explicitly broken.) Therefore, r2(A2, A3) is given to sextic order by 
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P2«2|A2p + PsaalAsP + - /32|A2|^ + /33I Agj^ + /332|A2nA- 



2|^ + 73|A3|^ + 7322|A3nA2|4 + ^^Zz]^^^]^^? 



(2) 



where we have chosen notation consistent with that used in the two flavor study of Ref. |H] , which arises as 
a special case of ^ if we take A2 or A3 to be zero. Pi is the number of vectors in the set {q^}. The form 
of the Ginzburg-Landau expansion ^ is model- independent, whereas the expressions for the coefficients a/, 
/3/Ji 7/1 and ^ijj for a specific crystal structure are model-dependent. For exchange symmetric crystal 
structures, 02 = as = /32 = /Ss = 72 = 73 = 7 and 7233 = 7322- 
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Because setting one of the A/ to zero reduces the problem to one with two-flavor pairing only, we can 
obtain a, (3 and 7 via applying the two-flavor analysis described in Ref. [S] to either {q2} or {q3} separately. 
Using a as an example, we learn that 

Here, qi = |q/| and A2SC is the gap parameter for the 2SC (2-flavor, 2-color) BCS pairing obtained with 
djjLi = and A/ nonzero with the other two gap parameters set to zero. Assuming that Aq <SC /x, it is given 
by A2SC = 23 Aq |1L)| . In the Ginzburg-Landau approximation, in which the A/ are assumed small, we must 
first minimize the quadratic contribution to the free energy, and only then investigate the quartic and sextic 
contributions. Minimizing a/ fixes the length of all the vectors in the set {q/}, yielding qi = "qS^i with 
rj = 1.1997 the solution to ^ log [(77 + !)/(?/ — 1)] = 1 (Sj. Upon setting qj = -q 5^i, ^ becomes 

Furthermore, the only dimensionful quantities on which the quartic and sextic coefficients can depend are 
then the (5/i/ [HIE], meaning that for exchange symmetric crystal structures and with 5^2 = (5/^3 = 5^ we have 
(3 = P/Sjj?, P32 = /332/'5/i^, 7 = 7/<5m^ and 7322 = 7322/(5 /i** where the barred quantities are dimensionless 
numbers which depend only on {€[2} and {€[3} that must be evaluated for each crystal structure. Doing 
so requires evaluating one loop Feynman diagrams with 4 or 6 insertions of A/'s. Each insertion of A/ 
(Aj) adds (subtracts) momentum 2qj for some a, meaning that the calculation consists of a bookkeeping 
task (determining which combinations of 4 or 6 q^'s are allowed) that grows rapidly in complexity with 
the compexity of the crystal structure and a loop integration that is nontrivial because the momentum in 
the propagator changes after each insertion. See Ref. 0], where this calculation is carried out explicitly for 
11 crystal structures in a mean-field NJL model upon making the weak coupling (Aq and Sfi both much 
less than fi) approximation. Note that in this approximation neither the NJL cutoff nor the NJL coupling 
constant appear in any quartic or higher Ginzburg-Landau coefficient, and they appear in a only within Aq. 
Hence, the details of the model do not matter as long as one thinks of Aq as a parameter, kept <^ /i. 

It is easy to show that for exchange symmetric crystal structures any extrema of Q{A2, A3) in (A2, A3)- 
space must either have A2 = A3 = A, or have one of A2 and A3 vanishing It is also possible to show that 
the three-flavor crystalline phases with A2 = A3 = A are electrically neutral whereas two-flavor solutions in 
which only one of the A's is nonzero are not 0j. We therefore analyze only solutions with A2 = A3 = A. In 
marked contrast to the two-flavor results of Ref. [H] which (ignoring the requirement of neutrality) show that 
many two-flavor crystal structures have negative 7 and hence sextic order free energies that are unbounded 
from below, we flnd that il{A, A) is positive for large A for all the crystal structures that we investigate 0]. 
This allows us to minimize ri(A, A) with respect to A, thus evaluating A and Q. 

We begin with the simplest three-flavor "crystal" structure in which {q2} and {q3} each contain only a 
single vector, making the (us) and (ud) condensates each a single plane wave. This simple condensate yields 
a qualitative lesson which proves helpful in winnowing the space of multiple plane wave crystal structures |1] . 
For this simple "crystal" structure, all the coeffflcients in the Ginzburg-Landau free energy can be evaluated 
analytically The terms that occur in the three-flavor case but not in the two-flavor case, namely (332 
and 7322, describe the interaction between the two condensates, and depend on the angle (j) between qj and 
q3. For any angle (j), both P32 a-nd 7322 are positive. And, both increase monotonically with and diverge 
as ^ TT 0]. This divergence tells us that choosing q2 and q3 precisely antiparallel exacts an infinite 
free energy price in the combined Ginzburg-Landau and weak-coupling limit in which A <^ Sfi, Aq <C /i, 
meaning that in this limit if we chose (/> = tt we find A = 0. Away from the Ginzburg-Landau limit, when 
the pairing rings on the Fermi surfaces widen into bands, choosing (j) ~ n exacts a finite price meaning 
that A is nonzero but smaller than that for any other choice of 0. The high cost of choosing q2 and q3 
precisely antiparallel can be understood qualitatively as arising from the fact that in this case the ring of 
states on the u-quark Fermi surface that "want to" pair with d-quarks coincides precisely with the ring 
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that "wants to" pair with s-quarks [5]. This simple two plane wave ansatz has been analyzed upon making 
the weak-coupling approximation but without making the Ginzburg-Landau approximation j^- AH the 
qualitative lessons learned from the Ginzburg-Landau approximation remain valid and we learn further that 
the Ginzburg-Landau approximation always underestimates A 

The analysis of the simple two plane wave "crystal" structure, together with the observation that in more 
complicated crystal structures with more than one vector in {q2} and {q3} the Ginzburg-Landau coefficient 
P32 (7322) is given in whole (in part) by a sum of many two plane wave contributions, yields one of two rules 
for constructing favorable crystal structures for three-flavor crystalline color superconductivity {q2} and 
{qg} should be rotated with respect to each other in a way that best keeps vectors in one set away from the 
antipodes of vectors in the other set. The second rule is that the sets {q2} and {q^} should each be chosen 
to yield crystal structures which, seen as separate two-flavor crystalline phases, are as favorable as possible. 
The 11 crystal structures analyzed in Ref. [5j allow one to make several pairwise comparisons that test these 
two rules. There are instances of two structures which differ only in the relative orientation of {q2} and {qg} 
and in these cases the structure in which vectors in {q2} get closer to the antipodes of vectors in {qg} are 
disfavored. And, there are instances where the smallest angle between a vector in {q2} and the antipodes of 
a vector in {qg} are the same for two different crystal structures, and in these cases the one with the more 
favorable two-flavor structure is more favorable. These considerations, together with explicit calculations, 
indicate that two structures, which we denote "2Cube45z" and "CubeX", are particularly favorable. 

In the 2Cube45z crystal, {q2} and {qg} each contain eight vectors pointing at the corners of a cube. 
If we orient {q2} so that its vectors point in the (±1,±1,±1) directions in momentum space, then {qg} is 
rotated relative to {q2} by 45° about the z-axis. In this crystal structure, the (ud) and {us) condensates 
are each given by the most favored two-flavor crystal structure The relative rotation maximizes the 
separation between any vector in {q2} and the nearest antipodes of a vector in {qg}. 

We arrive at the CubeX structure by reducing the number of vectors in {q2} and {qg}. This worsens 
the two- flavor free energy of each condensate separately, but allows vectors in {q2} to be kept farther away 
from the antipodes of vectors in {qg}. We have not analyzed all structures obtainable in this way, but we 
have found one and only one which has a condensation energy comparable to that of the 2Cube45z structure. 
In the CubeX structure, {q2} and {q^} each contain four vectors forming a rectangle. The eight vectors 
together point toward the corners of a cube. The 2 rectangles intersect to look like an "X" if viewed end-on. 
The color, flavor and position space dependence of the CubeX condensate is given by 



27r , , 2tt , . 

cos — [x + y + z) + cos — (— .T — y z) 
a a 



27r , , 27r , , 

cos — (— x + y + zj + cos — [X — y + z) 



(5) 



where a ~ ^fi-K/q = 4.536/(5/i = /.t/(1.764Mg ) is the lattice spacing. For example, with / i^i ~ 
100, 150, 200 MeV the lattice spacing is a = 72, 48, 36 fm. We depict this condensate in Fig. ^ 

In Figs. 2 and 3, we plot A and fl versus M^/^ for the most favorable crystal structures that we 
have found, namely the CubeX and 2Cube45z structures described above. We have taken the CFL gap 
parameter Aq = 25 MeV in these figures, but they can easily be rescaled to any value of Aq ^ ^ @]. 
Fig. 2 shows that the gap parameters are large enough that the Ginzburg-Landau approximation is at the 
edge of its domain of reliability. However, results obtained for simpler crystal structures suggest that the 
Ginzburg-Landau calculation underestimates A and the condensation energy and that, even when it breaks 
down, it is a good quahtative guide to the favorable structure (2]. We therefore trust the result, evident in 
Fig. 3, that these crystalline phases are both impressively robust, with one or other of them favored over a 
wide swath of M'^/fi and hence density. We do not trust the Ginzburg-Landau calculation to discriminate 
between these two structures, particularly given that although we have a qualitative understanding of why 
these two are favorable we have no qualitative argument for why one should be favored over the other. 
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Figure 1. The CubeX crystal structure 
of Eq. JFJ. The figure extends from to 
a/2 in the x, y and z directions. Both 
A2{r) and A3(r) vanish at the horizontal 
plane. A2(r) vanishes on the darker vertical 
planes, and A3(r) vanishes on the lighter 
vertical planes. On the upper (lower) dark 
cylinders and the lower (upper) two small 
corners of dark cylinders, A2(r) = +3. 3 A 
(A2(r) = — 3.3A). On the upper (lower) 
lighter cylinders and the lower (upper) two 
small corners of lighter cylinders, A3(r) = 
-3.3A (A3(r) = +3.3A). The largest value 
of I A/(r)| is 4A, occurring along lines at the 
centers of the cylinders. The lattice spacing 
is a when one takes into account the signs of 
the condensates; if one looks only at | A/(r)|, 
the lattice spacing is a/2. 
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We are confident that 2Cube45z is the most favorable structure obtained by rotating one cube relative to 
another. Wc are not as confident that CubcX is the best possible structure with fewer than 8+8 vectors. 
Regardless, the 2Cube45z and CubeX crystalline phases together make the case that three-flavor crystalline 
color superconducting phases are the ground state of cold quark matter over a wide range of densities. If 
even better crystal structures can be found, this will only further strengthen this case. 

Fig. 3 shows that over most of the range of M"^ / fjL where it was once considered a possibility, the gCFL 
phase can be replaced by a much more favorable three-flavor crystalline color superconducting phase. It 
also shows that it is hard to find a crystalline phase with lower free energy than the gCFL phase at the 
lowest values of Af^//i (highest densities) in the "gCFL window". This narrow window where the gCFL 
curve remains the lowest in Fig. 3 is thus the most likely place in the QCD phase diagram to find the gCFL 
phase augmented by current-carrying meson condensates described in Refs. Except within this window, 
crystalline color superconducting phases with the CubeX or the 2Cube45z crystal structures provide an 
attractive resolution to the instability of the gCFL phase. 

The three-fiavor crystalUine color superconducting phases with CubeX and 2Cube45z crystal structures 
are the lowest free energy phases that we know of, and hence candidates for the ground state of QCD, 
over a wide range of densities. One or other is favored over the CFL, gCFL and unpaired phases for 
2.9Ao < M'^/^i < 10.4Ao, as shown in Fig. 3. For Aq = 25 MeV and Ms = 250 MeV, this translates to 
240MeV < n < 847MeV. With these choices of parameters, then, the lower part of this range of n (higher 
part of the range of M"^ / fi) is certainly superseded by nuclear matter. And, the high end of this range 
extends far beyond the fi ^ 500 MeV characteristic of the quark matter at the densities expected at the 
center of compact stars. If compact stars have quark matter cores, then, it is reasonable to include the 
possibility that the entire quark matter core could be in a crystalline color superconducting phase on the 
menu of options that must ultimately be winnowed by confrontation with astrophysical observations. If Aq 
is larger, say ^ 100 MeV, the entire quark matter core could be in the CFL phase. 

Now that we have two candidates for the crystal structure of the three-flavor crystalline color 
superconducting phase of cold quark matter, favorable over a very wide range of intermediate densities, and 
a qualitative guide to the scale of A and fi(A), we can look ahead toward the calculation of astrophysically 
relevant observables. The heat capacity of the CubeX and 2Cube45z phases should be only quantitatively 
suppressed relative to that of unpaired quark matter, but their neutrino emissivity may be more significantly 
suppressed Both can now be calculated, yielding an estimate of the eff'ects of a crystalline quark matter 
core on the rate at which a neutron star cools by neutrino emission. A crystalline color superconducting 
core, being both crystalline and supcrfluid, could be a region within which rotational vortices are pinned and 
hence (some) pulsar glitches originate 0- Or, the presence of crystalline quark matter within neutron stars 
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Figure 2. Gap parameter A versus 
/ ^ for the CFL gap parameter, the 
three gap parameters Ai < A2 < 
A3 describing (ds), (us) and {ud) 
pairing in the gCFL phase, and the 
gap parameter in the crystaUino color 
superconducting phases with CubcX 
and 2Cube45z crystal structures. 



Figure 3. Free energy Q relative to 
that for neutral unpaired quark matter 
versus Mf/^, for the same phases as 
in Fig. 2. The gCFL phase free energy 
is a dashed line as a reminder that this 
phase is unstable: the free energy of the 
ground state of matter must lie below 
this dashed line. 



could be ruled out if it predicts glitch phenomenology in qualitative disagreement with that observed. The 
two microphysical properties of crystalline quark matter that must be estimated before glitch phenomenology 
can be addressed are the vortex pinning force and the shear modulus. Glitches occur if rotational vortices are 
pinned and immobile while the spinning pulsar's angular velocity slows over years, with the glitch triggered 
by the catastrophic unpinning of long-immobile vortices. Immobilizing vortices requires sufficient pinning 
force and shear modulus. Estimating the pinning force will require analyzing how the CubcX and 2Cube45z 
phases respond when rotated. The shear modulus is related to the coefficients in the effective theory that 
describes the phonon modes of the crystal Analyzing the phonons in the three-flavor crystalline color 
superconducting phases with the 2Cube45z and CubeX crystal structures is thus a priority. These are the 
prerequisites to determining whether observations of pulsar glitches can be used to rule out (or in) the 
presence of quark matter in a crystalline color superconducting phase within neutron stars. 

Wc thank the LBNL Nuclear Theory Group for hospitality. Research supported in part by the U.S. DOE 
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